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The universal density profile of the central region of dark 
matter haloes 
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ABSTRACT 

We consider the density profile of the central region of dark matter haloes. It turns 
out that under very general conditions the profile is universal: it depends almost not 
at all on the properties of the initial perturbation and is very akin, but not identical, 
to the Einasto profile. 

We estimate the size of the 'central core' of the distribution, i.e., the extent of the 
very central region with a respectively gentle profile, and show that the cusp formation 
is unlikely, even if the dark matter is cold. We also indicate that the density profile of 
the outer part (r > 0.5R V i r ) of the haloes significantly depends on the initial conditions 
and should not be universal, in contrast to the central area. All these results can be 
useful both to indirect search of the dark matter and to N-body simulations of the 
structure formation. 

Key words: cosmology: dark matter, cosmology: large-scale structure of Universe, 
elementary particles, methods: analytical. 



1 INTRODUCTION 

The problem of the Universe structure formation is still far 
from the complete solution. Since the task is essentially non- 
linear, N-body simulation is one of the possible ways to solve 
it. An exceptional result found in the simulations is that ex- 
tremely varied initial density profiles evolve very soon into 
a very similar 'universal' density profile (hereafter UDP). 
The UDP closely resembles the Na varro-Frenk-White pro- 
file (|Navarro. Frenk. fc White! Il997l), though further simu- 
lations show llMoore. Lake, fc Katdll998l: fcao et al.11200 ' 
iDiemand et al.ll2008l ; IStadel et al.ll2009l:lNavarro et al.ll2oT 



that it is closer to the Einasto one l|Einastdll969E i.e., the 
density in the very centre of the halo tends to a constant 
value instead of behave as r^ 1 . 

The reasons of the universal density profile appearance 
are still unclear. It seems likely that the presence or absence 
of substructures is not determinative for the process: the 
UDP appears, even if we consider a halo col lapse with no 
substructures (iHuss. Jain, fc Steinmet j fl999h or hot dark 
matter haloes jWang fc Whitd l2009fT of course, N-body 
simulations are not infallible: in particular, they model mi- 
nuscule dark matter particles by 10 4 — 10 6 A/q bodies. Con- 
trary to the particles, the bodies can scatter on large angles 
during close approaches, and in order to avoid this effect, 
the newtonian potential of test bodies is modified. However, 
the N-body simulations excellently describe observed prop- 
erties of galaxy clusters; in particular, their density profiles 



are close to the UDP. Therefore the modelling results are ad- 
equate, at least, in the important case of the highest mass 
systems, which indicates that the universal density profile is 
not a calculational artefact and has s ome physical motiva- 
tion l|Dalal. Lithwick. fc KuhlerJl201Ch . 

It is interesting to note, however, that if we consider 
a scale-free initial profile, UDP often doe s not appear, and 
the fi nal profile depend on the initial one l|Lithwick fc Dalall 
l201ll) . We should mention that such profiles cannot describe 
the structure formation in the real Universe, where there is 
always some characteristic initial radius of the region, 
from which the majority of the particles will be captured 
by the forming halo, while the matter outside of JH will be 
captured by the other nearby structures or form a diffuse 
component. As we will see below, the existence of charac- 
teristic scale of the initial profile plays an important role 
in the UDP formation. 

We do not try to fully explain the origin of the uni- 
versal density profile in this letter. We shall show, however, 
that under some very general conditions, the central part of 
the formed halo should really have almost the same density 
profile, close, but not equal, to the Einasto one. 



2 ENERGY SPECTRUM 

In order to explain some properties of the distribution 
of real systems, the 'violent relaxation' is frequently in- 
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volved, i.e. the intensive energy exchange between dark 
matter particles at the initial stage of the collapse, when 
strong small-scal e gravitational fie lds occur in the centre 
part of the halo |Lvnden-BeUlll967l ). Further investigations 
showed, however, that the situation is more likely quite the 
contrary: the energy of the majorit y of the particles re- 
mains quite close to the initial one ijOuinn fc Zureklll98Sl ; 
IZaroubi. Nairn, fe Hoffmanll 19961 ). It is the fact that we take 
as the starting point: it means that the final energy distri- 
bution of the particles is respectively close to the initial one. 

Let us consider the initial distribution of matter. Here- 
after in this article we will consider a spherically-symmetric 
task and neglect the presence of baryon matter. The mass 
and initial radius of the halo are m and JR. Let us denote 

2 

by e = + <f> the total energy of a unit mass of the matter, 
(f> is the gravitational potential. If the dark matter is cold, 
the initial velocity of the matt er may be thought of as to b e 
zero without loss of generality l)Gorbunov fe Rubakovl201ll ) . 
(As we will see at the end, the conclusions of this section is 
also valid for the warm dark matter). To begin with, let us 
consider the case when the matter is uniformly distributed 
inside 5H (p = const, the Tolman case). Inside the sphere: 

♦m-SG-Kh) 1 ) « 

Hence we obtain the initial energy distribution of the sub- 
stance: 

where $ = ^M-- Performing similar calculations, we obtain 
for the case of initial density perturbation par -1 

- = -(2 + ¥ ); eg [-»;-*] (4) 

and for initial profile p oc r~ 2 

dm m ( e \ . . . 

$ exp ( 1+ $) ; ^[-oo;-*] (5) 

As one can see in Fig.[T] distributions ©, Q, and even ^ 
differ not so strongly. 

It seems reasonable to say that the initial energy dis- 
tribution is similar to Q under any reasonable choice of 
shape of the initial perturbations. Indeed, a real (smooth) 
density perturbation collapses when — = /? ~ 1, and /3 
weakly depends on the shape of the initial perturbations. 
Consequently, if the initial shape of 8p(r) is a smooth func- 
tion, distributions of <f> and dm differs from {TJ and ((2]) only 
on a multiplier of the order of 2. Consequently, distribution 
dm/de does not differ very significantly from the uniform 
one ([3]). Second, if we divide a spherically symmetric initial 
distribution on radial layers of thickness dr, e is constant 
within the layer, and the volume of the layer grows as (=j) . 
As a result, the particles with e ~ — $ (which corresponds to 
r = 9V) dominate in the initial energy spectrum even in the 
case of quite a steep profile (8p ~ r _1 ). A strong departure 
from (O will take place only if the initial profile is very steep 
(for instance, Sp ~ r~ 3 ). On the other hand, such profiles 
are hardly suitable to describe the initial perturbations of 
the dark matter. 
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Figure 1. Energy spectra for various forms of initial pertur- 
bations. Apparently, there are no particles with e > — 



So for a large variety of initial conditions the initial 
energy distribution of the particles is similar to (|3]), at least, 
in the following sense: the majority of the particles have the 
energy close to — $, and their energy distribution is quite 
narrow. In case of distribution (|3}, for instance, for all the 
particles e £ [— | $;—$]. This conclusion is also valid for 
the warm dark matter. Indeed, the smearing of the energy 
spectrum with respect to the cold dark matter case is of 
the order of the temperature of dark matter at the moment 
of the structure collapse. For the structure to be able to 
collapse, the temperature should be much lower than $. 

As we have already mentioned, the energy distribution 
of the particles changes quite weakly with respect to the ini- 
tial one. Therefore, the similarity between the initial distri- 
butions leads to the similarity of the particle energy spectra 
of the relaxed halo. But how can the system evolve at all, if 
the energies of the bulk of the particles remain almost the 
same? The mechanism is the following: the particles oscil- 
late in the potential well. Initially they oscillate cophasal. 
Then, even if the energies of the particles remain constant, 
their phases diverge. Typically, the period sharply depends 
on e for the majority of potential wells. But even if the en- 
ergies of the particles precisely coincide, their phases can di- 
verge as a result of a small nonsphericalness, perturbations 
from the nearby structures etc. Finally, the phases of the 
oscillations are mixed thoroughly, and each phase becomes 
equiprobable. This is exactly the situation, corresponding 
to the stationary halo. It is important to underline that the 
phase mixing goes on continuously, in contrast to the energy 
exchange, which is effective only during the first short stage 
of violent relaxation. 

It is significant that the small energy evolution of the 
majority of the particles is not only an effect discovered by 
N-body simulations, but it has deep physical reasons. As we 
could see, the main fraction of mass of the future halo is 
initially concentrated in the outer layers with the highest 
initial r (since thei r volume grows as (r/9l) 2 ). Meanwhile, 
as it was shown by iLvnden-Belll (|l967l ). the efficiency of all 
relaxation processes (and so the alternation of the particle 
energy) rapidly drops with r: even in that work the outer 
part of the system remained unrelaxed. 
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However, the strongest argument in favour of the energy 
spectrum conservation during the relaxation is the existence 
of the standard density profile itself. Indeed, let us suppose 
that a strong energy trade between dark matter particles 
occurs. Then the energy exchange can appear only on the 
initial stage of the halo formation: the energy exchange is 
implemented by the variable small-scale gravitational fields 
appearing only just after the halo collapse. It is apparent, 
however, that in this case the energy trade strongly depends 
on the initial conditions. Indeed, if the initial perturbation 
is flat (above-mentioned Tolman case ©), all the particles 
fall on the centre simultaneously. On the contrary, in the 
case of p ~ r _1 or p ~ r~ 2 initial profiles the outer particles 
fall on the centre much later, than the inner ones. More- 
over, as long as the outer particles fall to the centre, the 
inner ones have already oscillated many times. So when the 
outer dark matter passes through the central region the first 
time, the region has already relaxed. It is clear that the con- 
ditions of the energy trade are completely different in this 
two cases, and if the density profile is defined by the energy 
redistribution, we might expect that the final density pro- 
files were completely different. Since even though the final 
density profiles are similar, it means that the influence of 
the energy redistribution is not principle. 



3 DENSITY PROFILE 

Finally, the halo collapse forms stationary potential well 
4>(r). Each particle can be characterized by maximum ra- 
dius ro the particle moves from the centre. As we could see 
©, the initial energies of the particles are close, and the 
energy distribution does not change much during the relax- 
ation. Consequently, ro of the majority of the particles are 
concentrated in a fairly narrow interval below some charac- 
teristic radius R (R is comparable, but does not obligatory 
coincide with fR). The fractions of the halo mass formed by 
particles with ro > R or ro <C R are small. It follows that 

R — Rvir - 

On the other hand, it means that the main contribution 
to the density in the central part of the halo is made by the 
particles that come there from the outside. Indeed, a formed 
halo has 

k = \{ct>{R vlr ) - <t>(0))/<t>(Rvir) - <t>(oo)))\ > 1 

For instance, the Navarro-Frenk- White halo h as k = 
Cyjr / lnfcw + !) — !, i.e. k — 1.5 even if cw = 3 jBaushev! 



3,0 



|2012| ). Comparing this with distribution we may con- 
clude that the majority of the particles has ro > 0.5R V i r - 
Of course, the real distribution differs from and there 
are always some particles, which energies have completely 
changed, and their ro < 0.5R V i r - However, more than the 
half of the halo mass is typically concentrated inside r = 
0.5R V ir, and if the particles with ro < 0.5R V i r dominated 
there, it would mean that the energies of approximately the 
half of the halo particles have changed several times with 
respect to the initial ones, which is highly improbable, as 
we could see. 

We can easily find the density distribution in the cen- 
tre of the halo, created by the particles coming to this re- 
gion from the outside. Let us single out the particles with a 
certain ro. We denote their total mass by m. According to 
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Figure 2. Ratios T(/x)/T(0) for the newtonian potential <f> oc 
—X/r (solid line), isothermal sphere <f> oc ln(r) (dashed line) and 
Navarro-Frenk-White profile with c v i r = 5 (dotted line). The 
normalization of the potentials is chosen so that the mass inside 
ro is the same for all three cases. p c i rc corresponds to the circular 
orbit. 



results of the numerical simulations, we assume that their 
specific angular momentum p = [v x r\ has Gaussian distri- 
bution 



2p 



dm — m^r exp ^ du 

~v2 \ n/^ ' 



(6) 



A particle moving in gravitational field c/>(r) has two inte- 

2 2 
V Li 

grals of motion: p and e = — + ^-^ + tb. The radial velocity 
of the particle is equal to 



tv = W2Mro)-#r))-Ai 3 (3-i 



(7) 



We introduce the maximum angular momentum of a particle 
wherewith it can reach radius r 



/ 1 

Umax = 2(0(r O ) ~ Cf>(r)) ( ^ 2 



Then (Q) may be rewritten as 



rro 



I'- 



(8) 



(9) 



A particle having maximal radius ro and minimal radius 
r-min contributes into the halo density on all the interval 
[r m in, ro]. The contribution of a single particle of mass m on 
an interval dr is proportional to the time the particle passes 



an interval ar is proportional to t 
on this interval jBaushevll20Tll ). 



dm 
m 



dt 
T 



dr 



(10) 



Here T is the half-period of the particle, i.e. the time it 
takes for the particle to fall from its maximal radius to the 
minimal one 

T(r„,n)=/ - (11) 

Jr min IV 

T(ro, At) is, generally speaking, a function of ro and fi. 
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Figure 3. The density profile of the exact solution of (1 1 3 I t for 
f(ro) = M v i r S(ro — R v ir) (dashed line) and profile l|19j l with 
the same r c /R v i r (solid line). One can see that the departure 
of approximative equation (| 1 9 1> from the exact solution is quite 
small. 



We can obtain the density distribution from (|10[1 . sub- 
stituting there ((9} instead of v r and 4nr 2 p instead of dm: 



mro 



4-Kr 2 v r T AnrTs/rl - r 2 y/p 



(12) 



In order to find the density distribution produced by all the 
particles of the halo with a certain ro (their total mass we 
have denoted by m), we should integrate (|12[) over distribu- 
tion (©. 



P = 



mro 



2/zexp (— p 2 /a 2 ) 



Auryjr 2 -r 2 J a 2 T(r ,p)^p 2 nax - p} 



dp (13) 



This is the exact solution; we can significantly simplify 
it, however, if we take into account that T(ro,p) is in gen- 
eral a very weak function of p, especially for small p, since 
T(ro, p) has an extremum at p = 0. As we see in Fig.(2] T(p) 
differs from T(0) on no more, than 25%, for any reasonable 
potential, if p ^ 0.5p C irc, where p C irc is the angular momen- 
tum corresponding to the circular orbit. Only particles with 
small p can reach the central region under consideration: 
consequently, we may approximate T by T(ro) = T(ro,0): 
Then we can take the integral in (| 13p 



mro 



2naT(ro)ry / rg — r 2 



D 



ftSSS.^ (14) 



where D is the Dawson function D{x) = e~ x J Q X e* dt. Of 
course, all particles in a real halo have some distribution 
over ro: 



dm = f(r )dr 



(15) 



In order to take it into account, we should substitute 
/(ro)dro instead of m to (|14p and integrate the result over 
ro- Moreover, equation (|14[) can be significantly simplified 
in the central part of the halo. First of all, gravitational 
potential of the centre of the halo, 0(0), is finite for any 
reasonable halo profile (otherwise the annihilation signal, 
produced by the halo, would be infinite, being proportional 



to J 4-ivr 2 p 2 dr). Then we may rewrite iJSJ as 

2 n l(M \ J.tn\\ (i - 0(0) ^ rg 

p max = 2r (0(r o ) - 0(0)) [l - ^ _ m ) 

(16) 

The two last factors may be neglected in the center 
of the halo, and we can use approximations p ma x — 
r^/2((j>(ro) — 0(0)), \Jr\ — r 2 ~ ro. So the total density dis- 
tribution in the central part is 



POO 

Jo 2na 



/(ro) 



(r )T(r )r 



1 r ( 7) 



Here we use designation a(ro) instead of a in order to em- 
phasize that a is, generally speaking, a function of ro. 

As we could see, function /(ro) sharply depends on 
ro and differs noticeably from zero only on the inter- 
val between ~ Rvir and ~ R V i r /2. On the contrary, 
a(ro) should not change much on interval [Rvir, Rvir/ 2]: 
a(ro) is widely believed to be a powe r-law depen- 
dence with the index between —1 and 1 |Hansen et al.l 
2006). ^2(0(7-0) -0(0)) chan ges even slower: for instance, 
vWW) " 0(O))/(0(i?„ lr /2) - 0(0)) ~ 1.13 for the NFW 
profile with c„ ir = 5. Comparing this with the sharp be- 
havior of /(ro), we may neglect the weak dependence of the 
argument of function D in (|17[) on ro and substitute there 
some value, averaged over the halo 



a(r ) 



*(r )) 



v /2(0(r o )-0(O))/ 
Then we can rewrite (|17p and obtain the final result: 

1 f°° f(ro)dr 



rc^fr 
P = Pc— D[ — 

r \r c 



Pc = 



2%r c 



(18) 



(19) 



/0 Q(ro)T(r ) 

Since D(r/r c ) ~ r/r c , when r/r c — > 0, p c is the central 
density of the halo. As we can see, it is always finite. At the 
same time, the shape of the density profile depends on the 
only parameter r c . 



4 DISCUSSION 

Before proceeding any further, we should estimate the char- 
acteristic value of a. This question is, generally speaking, 
not quite clear. According to the modern cosmological con- 
ceptions, the initial ang ular momentum of each dark matter 
particle was very small (|Gorbunov fc Rubakovll201lh , and it 
is the assumption that is routinely used in the N-body simu- 
lations. However, the particles can acquire some momentum 
in the subsequent relaxation. There are some grounds for 
believing that the mean-square-root angular momentum of 
the dark matter in the h alo of the Milky Way, for instance, 
is small (|Baushev|[20Tlh . A detailed consideration of this 
problem goes far beyond the scope of this letter. However, 
N-body simulations generally result in the mean-square-root 
angular momentum, close to the maximum possible one, 
and, as our consideration is oriented to the numerical calcu- 
lations, we will proceed from this assumption and accept 



(a) 



- Rvir 



Rvir 



— yj G Ad^v i r Rv i r 

4 



(20) 



which roughly corresponds to the situation, when all par- 
ticles with ro = 0.75iiW and v < v eac at this radius fall 
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Figure 4. Density profile H19H of the central region of a halo with 
r c = _R„i r /5.74. Einasto profile with n = 0.5 and r s = 0.06R v i r 
is plotted for comparison (dashed line). 

within the 3cr interval of distribution (|6j). It is apparent that 
(a) cannot substantially exceed (|20|) : otherwise a significant 
part of the dark matter would not be gravitationally bound 
in the halo. 

In deriving (|19|l we used several times the fact that r <C 
ro. This raises the question: for what r/R V i r is equation 
(|19[) still valid? In order to check this, we solved numerically 
exact equation (|13[) on condition ((6J and for distribution 
/(ro) = M vir 8(ro — R V i r ). 4>(r) appears in expression (fl3)) . 
and we need to write out equation 

d<t> g r. 2 , s 

— = — Atvx p(x)dx 
dr r 2 J 

in order to close the task. The exact density profile (dashed 
line) and profile (|19|) with the same r c /R V i r (solid line) are 
represented in Fig. [3] We see that (|19p well approximates 
the exact solution up to r ~ 0.3R V i r - So (|19p can be valid 
in quite a large area around the halo centre. 

Fig.[4]represents density profile (|19|l for r c = _R„i r /5.74. 
For reference, Einasto profile with n — 0.5 and r s = 0.06R V i r 
is also plotted. As we can see, the profiles are as a whole 
similar. However, (|19[) is steeper than the Einasto profile, 
and in this sense is closer to the Navarro-Frenk- White one. 
Taking into acco unt large uncertai nties of the N-body cal- 
culation results l| Stadel et alj|2009h . we may conclude that 
profile (|19p is in excellent agreement with the simulations. 

In conclusion, we would like to summarize the main con- 
clusions of this letter. First, the density profile in the central 
part of a formed halo is really quite universal and weakly 
depends on the shape of the initial perturbation. This is a 
result of the fact that the potential well of the formed halo 
is much deeper than that of the initial perturbation. As a 
result, quite a large area builds up in the centre of the halo, 
where the particles do not reside, but just come there from 
the outside: their apoapsides lie outside of the area. In order 
that the orbit of a particle entirely lie within the area, the 
particle energy should be dramatically decreased with re- 
spect to the initial value. Such a process is always hampered 
in a dissipationless system, and hence the particles, coming 
to the central area from the outside, dominate there and 



form Einasto- like universal density profile (|19[) . The shape 
of the profile is totally determined by single parameter r c 

Second, the density in the centre of a halo should be 
finite even in the case, when the dark matter is cold. In- 
deed, one can see that equation (|13[l always asymptotically 
converges to (|19p when r — > (i.e., the central density re- 
mains finite), even if the angular momentum distribution 
completely differs from (|6} . There are only two possible sit- 
uations, when this conclusion is incorrect, and the cusp for- 
mation is possible. It is feasible if angular momentum distri- 
bution © has a peculiarity at [i = 0, i.e. the halo contains an 
anomalous number of particles with virtually zero angular 
momentum. This seems very unlikely. The cusp formation 
is also possible, if the halo contains a lot of particles with 
the energies, almost precisely equal to </>(0). The particles 
continuously reside in the centre of the halo and it is these 
particles that form the cusp. It is worthy of note, however, 
that the mass of the particles should be quite significant. In 
the case of Navarro-Frenk- White profile with cw = 5—10, 
for instance, the cusp area (with the profile, close to r -1 ) 
comprises 13 — 20% of the total halo mass. As it has been 
already mentioned, the energy of the particles should be 
at least several times decreased with respect to the initial 
value, and so wide variations of such a large fraction of the 
halo particles are improbable. Moreover, as we discussed in 
the energy spectrum section, the parameters of the energy 
exchange process critically depend on the initial conditions, 
i.e. the form and even the existence of the cusp should be 
defined by the shape of the initial perturbation. N-body sim- 
ulations indicate the opposite. 

Third, contrary to the central part, the density profile 
of the halo outer regions essentially depends on the shape 
of distribution /(ro) (|15p . /(ro) is, in its turn, sensitive to 
the details of the relaxation process. Consequently, we may 
expect that the density profile of the outer regions of haloes 
(r > 0.5R V ir) is not universal and depends on the initial 
conditions to a degree. 
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